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We consider a branching Brownian motion with linear drift 
in which particles are killed on exiting the interval (0,K) and 
study the evolution of the process on the event of survival as 
the width of the interval shrinks to the critical value at which 
survival is no longer possible. We combine spine techniques and 
a backbone decomposition to obtain exact asymptotics for the 
near-critical survival probability. This allows us to deduce the 
existence of a quasi-stationary limit result for the process condi- 
tioned on survival which reveals that the backbone thins down 
to a spine as we approach criticality. 

This paper is motivated by recent work on survival of near crit- 
ical branching Brownian motion with absorption at the origin 
by Ai'dekon and Harris in [AH] as well as the work of Bcrestycki 
et al. in [BBS] and [BBS11]. 

1. Introduction and main results. 

1.1. Introduction and main results. We consider a branching Brow- 
nian motion in which each particle performs a Brownian motion with 
drift —fi, for \i > 0, and is killed on hitting or K. All living parti- 
cles undergo branching at constant rate /3 to be replaced by a random 
number of offspring particles, A, where A is an independent random 
variable with distribution {q^ k = 0, 1, ...}, finite mean m > 1 and such 
that E(Alog + A) < oo. Once born, offspring particles move off inde- 
pendently from their birth position, repeating the stochastic behaviour 
of their parent. 

In other words, the motion of a single particle is governed by the in- 
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finitesimal generator 

defined for all functions u G C 2 (0,K), the space of twice continu- 
ously different iable functions on (0, K), with u(0+) = u(K-) = 0. The 
branching activity is characterised by the branching mechanism 

(1.2) F{s)=P(G(s)-s), sg[0,1], 

where G(s) = YlkLo 1kS k is the probability generating function of A. 
Let us introduce some notation. Denote by N t and |A^| the set of and 
the number of particles alive at time t respectively. For a particle u G 
Nt, we write x u (t) for its spatial position at time t. We define X t = 
YlueN ^x u (t) to be the spatial configuration of particles alive at time t 
and we set X = (X t , t > 0). Denote by P^ the law of X with X = v 
where v G A4 a (0,K), the space of finite atomic measures on (0,K) 
of the form ^" =1 S Xi with Xi G (0, K) and n G N. If the process is 
initiated from a single particle at x G (0,K), then we simply write 
(instead of P£)- We will sometimes neglect the dependence on the 
initial configuration and write P K without a subscript. We call the 
process X a P x -branching diffusion. 

Further, (£ = (£ t , t > 0), ) will henceforth denote a Brownian motion 
with drift — /i starting from x G (0, K) which is killed upon exiting the 
interval (0, K). F K is the law of the single particle motion under P K . 
For x G [0, K] we define the survival probability Pk( x ) — PxiC = °°) 
where ( = inf{t > : \N t \ = 0} is the time of extinction. As a first 
result we identify the critical width K below which survival is no longer 
possible. 

THEOREM 1. If ^ < a/2(?ti - and K > K where K : = 
7r(^/2(m — — yU 2 )" 1 , then Pk(x) > for all x G (0,K); otherwise 
p K (x) = for all x G [0, K]. 

The proof of Theorem 1 uses a spine argument, decomposing X into 
a Brownian motion conditioned to stay in (0, K) dressed with indepen- 
dent copies of (X, P K ) which immigrate along its path. 
As we want to study the evolution of the P x -branching diffusion on 
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survival, we will develop a decomposition which identifies the particles 
with infinite genealogical lines of descent, that is, particles which pro- 
duce a family of descendants which survives forever. To illustrate this, 
in a realisation of X, let us colour blue all particles with an infinite line 
of descent and colour red all remaining particles. Thus, on the event of 
survival, the resulting picture consists of a blue tree 'dressed' with red 
trees whereas, on the event of extinction, we see a red tree only. 
The branching rates of the branching diffusions corresponding to the 
blue and red trees can be intuitively derived as follows. For simplicity, 
consider the dyadic branching case only. A particles dies and is replaced 
by two offspring, at position x say, at rate (3. The probability that one 
of its offspring has an infinite genealogical line of descent is pk{x), in- 
dependent of the other particle. Thus, with probability Pk{x) 2 both off- 
spring particles are blue and likewise with probability (1 — Pk(%)) 2 and 
probability 2pk(x){1 —pk(x)) two red ones, respectively one blue and 
one red particle are born. Thus, given a particle is blue, it branches into 
two blue particles at rate /3 P pK\ = Ppk(-) and, given a particle is red, 
it branches into two red particles at rate ,5(1— Pk{')) while, given a par- 
ticle is blue, immigration of a red particle occurs at rate 2/3(1 — p K (-)). 
Similar reasoning gives the result for the general branching mechanism 
case, see Section 3. There we shall also see that particles in branch- 
ing diffusion corresponding to the red trees, respectively the blue tree, 
move according to a Brownian motion with drift — (/i + 1 ^ K K ) , respec- 
tively — (/!— ^). This results from /i-transforms of L using h — 1 —px 
and h = px respectively. In fact we will show that the laws of the 
branching diffusions corresponding to the red tree, the blue and the 
dressed blue tree arise from martingale changes of measure which, on 
the level of infinitesimal generators, correspond to the aforementioned 
^-transforms. 

Suppose we know the branching mechanisms of the branching diffusions 
corresponding to the blue and the red trees in the general well 
as the immigration rates (we will see that a branching mechanism of 
the general form induces a second type of immigration at the branching 
times of the blue tree). Intuitively speaking, the coloured tree starting 
from x G (0, K) is then constructed by flipping a coin with probability 
1 — pk{x) of 'heads' and if it lands 'heads' we grow a red tree with 
initial particle at x, while if it lands 'tails' we grow a blue tree at x and 
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dress its branches with red trees. Let us write for the law of the 
coloured tree which is defined on the filtration 3F% := o{T t -, c(u) u£ N t }, 
where (J-" t; t > 0) is the natural filtration of (X, P K ) and c(u) is the 
colour of particle u E N t . Note that the colours of all particles are F^- 
measurable. Then we can state (what will turn out to be a simplified 
version of) the so-called backbone decomposition. 

THEOREM 2 (Backbone decomposition). Let K > K and x E 
(0,K). On the filtration (.F^t > 0) (which means we ignore the colour- 
ing), (X, P^) is equal in law to (X, P^). This is, for all t £ [0, oo] and 
A E T t , we have Pf (A) = P?{A). 

An alternative statement of this result is given as Theorem 6 in 
Section 3. There we show that the backbone decomposition arises nat- 
urally from combining changes of measure which condition (X, P K ) on 
either the event of survival or the event of extinction. A significant con- 
venience of the backbone decomposition is that conditioning the P K - 
branching diffusion on survival is the same as conditioning on there be- 
ing a dressed blue tree, that is a blue tree 'dressed' with red trees. Thus, 
instead of studying the quasi-stationary limit lim^Xo Pxi'\C = °°) ^ 
suffices to study the evolution of the branching diffusion corresponding 
to a dressed blue tree as K I K . 

In order to do this, we need to know the asymptotics of the survival 
probability px near criticality. For a first asymptotic result note that 
u = 1—px solves the differential equation Lu + F(u) = on (0, K) with 
boundary condition u(0) = u(K) = 1 (cf. Remark 2). Near criticality 
we may assume that Pk{x) is very small for a fixed x and neglect- 
ing all terms of order {pk{x)) 2 and higher we obtain the linearisation 
LpK + mflpK = 0. This suggests Pk{x) ~ Ck sm(irx / K Q )e^ x . In fact 
we have the following result. 

Theorem 3. Uniformly for all x E (0, K ), we have 
(1.3) p K (x) ~ C K sm(<7rx/K )e» x , as K I K , 

where Ck is independent of x and can explicitly be determined as 

and in particular Ck 4 as K \, K . 
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In Section 4.1 we will prove a first part of Theorem 3, that is equa- 
tion (1.3), in the fashion of [AH] using spine techniques. It is to be 
particularly emphasized that we are able to determine Ck here. Ck is 
a 'non-linear' constant and 'linear' spine techniques fail when trying 
to identify it. Nevertheless, a careful application of the backbone de- 
composition given in Theorem 2 will deliver an explicit expression here 
since the blue tree captures enough 'non-linear' branching information 
about the evolution of (A, P ) on survival. A heuristic argument is 
given in Section 4.2.1 followed by the rigorous proof in Section 4.2.2. 
With Theorem 2 and 3 in hand we look for a quasi-stationary limit 
result for the law of the branching diffusion corresponding to a dressed 
blue tree, which agrees with the law of A conditioned on survival, as 
we approach criticality. In the dyadic case, the heuristic derivation of 
the branching rates already suggests that, given the particle positions 
x u (t) for u £ Nt in Xt, the number of particles in the blue tree at time 
t, is the number of successes in a sequence of independent Bernoulli 
trials each with probability of success pK(x u (t)), u G N t (We will ad- 
dress this thinning argument rigorously in Remark 1). Now, as K \. K , 
the probabilities Pk{-) tend to uniformly by Theorem 3 and thus the 
blue tree becomes increasingly thinner on (0, A ). Under conditioning 
on survival, it cannot vanish completely though since the genealogical 
line of the initial blue particle cannot become extinct and thus one may 
believe that, over a fixed time interval [0,T], the blue tree thins down 
to a single genealogical line at criticality. In the case of a dyadic branch- 
ing mechanism this conjecture can readily by confirmed by looking at 
the branching rates. The blue branching rate (3px drops down to as 
K 4- K , at the same time the red branching rate /3(1 — px) increases 
to (3 and the rate of immigration 2/3(1 — pk) rises to 2/3 at criticality. 
Formalising this idea and taking into account the change in the single 
particle motion, the general results reads as follows. 

Theorem 4. Let x e (0, K ). Consider a process X* = (X t *, t > 0) 
which evolves as follows. X* is initiated from a single particle at x 
performing a Brownian motion conditioned to stay in (0,K ), i.e. a 
strong Markov process with infinitesimal generator 
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defined for all u G C 2 (0, Kq). Along its path we immigrate A indepen- 
dent copies of (X, P K ) at rate m(3 where A has the size-biased offspring 
distribution (q^, k = 0, 1, ...) with 

k + 1 

q k = q k+ i , k > 0. 

m 

Denote the law of X* by P*. Then, for any fixed time T > 0, the law 
of (Xt, < t < T) under the measure lim^-^ P^{-\C = °o) is equal to 
(X*,0<t< T) under P*. 

The construction of the backbone via a martingale change of mea- 
sure allows us to give a very simple proof of this quasi-stationary limit 
result. Theorem 4 can be seen as an extension of the spine decomposi- 
tion we mentioned in the discussion following Theorem 1 to the critical 
width K . We emphasize however that the result, as stated, only holds 
over finite time horizons [0,T]. 

In Section 6, we demonstrate the robustness of our approach by apply- 
ing the results for the P^-branching diffusion to study the evolution of 
a supercritical super-Brownian motion with absorption at and K near 
criticality. We outline a backbone decomposition analogous to Theo- 
rem 6 in which we will see that the backbone of the super-Brownian 
motion with absorption at and K is the same as the backbone of an 
associated P^-branching diffusion. This connection allows us to deduce 
asymptotic results for the survival rate of the super-Brownian motion 
with absorption on (0, K) directly from the results on the survival prob- 
ability of the associated P x -branching diffusion. Further, we can find a 
quasi- stationary limit result for the super-Brownian motion equivalent 
to Theorem 4. This section is intended to highlight the applicability of 
the backbone approach and we will only sketch the proofs of the results 
therein. 

Our paper is organised as follows. In Section 2 we present the proof 
of Theorem 1 using spine techniques. In Section 3 we show that the 
backbone arises from a martingale change of measure which conditions 
(X, Pjf ) on survival, and we establish the backbone decomposition. We 
prove the asymptotic results for the survival probability given in The- 
orem 3 in Section 4. The proof of the quasi-stationary limit result in 
Theorem 4 follows in Section 5. Section 6 sketches the analogous results 
for the super-Brownian motion on (0,K). 
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1.2. Literature overview. Spine techniques of the type used in the 
proof of Theorem 1 were developed in Chauvin and Rouault [CR88], 
Lyons [Lyo97] and Lyons et al. [LPP95] and are now a standard ap- 
proach in the theory of branching processes. See, for example, Harris 
et al. [HHK06] and Kyprianou [Kyp04] for related applications in the 
setting of branching Brownian motion with absorption at respectively 
absorption at a space-time barrier. 

A backbone decomposition as in Theorem 2 for supercritical super- 
processes is presented in Berestycki et al. [BKMS11]. It extends the 
earlier work of Evans and O'Connell [E094], Fleischmann and Swart 
[FS04] and Englander and Pinsky [EP06] as well as the corresponding 
decomposition for continuous-state branching processes in Duquesne 
and Winkel [DW07]. 

The results for sup erpro cesses are complemented by the decomposi- 
tion in Etheridge and Williams [EW03] which considers the (1 + (3)- 
superprocess conditioned on survival. This work is of particular interest 
in the current context since it also presents the equivalent result for the 
approximating branching particle system. However we should point out 
that in their case the immigrants are conditioned to become extinct up 
to a fixed time T whereas, in our setting, we condition on extinction 
in the strip (0, K). Thus the underlying transformations in [EW03] are 
time-dependent in contrast to the space-dependent /i-transforms we see 
in our setting. 

We also point out that our derivation of the backbone decomposition 
differs from the previously mentioned articles in that we show that the 
backbone arises from combining changes of measure which condition 
(X, Pjf- ) on either the event of survival or the event of extinction. 

The equivalent result to Theorem 3 in the setting of Branching Brow- 
nian motion with absorption at the origin was shown in Berestycki et al. 
[BBS 11] and Aidekon and Harris [AH]. However, it has not been possi- 
ble so far to give such an explicit expression for the constant analogous 
to C K . 

A similarly fashioned result to Theorem 4 was obtained in the afore- 
mentioned work by Etheridge and Williams [EW03]. Their result ex- 
tends the Evans immortal particle representation for sup erpro cesses in 
[Eva93] which is the equivalent of the spine representation for branching 
processes. Again we point out that, in contrast to our setting, extinc- 
tion is a time-dependent phenomenon. Further, our martingale change 
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of measure approach to the backbone decomposition allows us to give 
a very simple proof of the quasi-stationary limit result. 

2. Spine techniques - Proof of Theorem 1. The proof of Theo- 
rem 1 uses classical spine techniques developed in Chauvin and Rouault 
[CR88], Lyons et al. [LPP95] and Lyons [Lyo97]; see e.g. Harris et al. 
[HHK06] and Kyprianou [Kyp04] for related applications in the setting 
of Branching Brownian motion with absorption at 0. 
We will briefly recall the key steps in the spine construction. For a com- 
prehensive account we refer the reader to Hardy and Harris [HH09]. 
Recall that we denote by (£, P£) a Brownian motion with drift —fj, 
initiated from x G (0, K) which is killed upon exiting (0, K). Then the 
process 

(2.1) T K (t) = sm{-Ki t /K)e^ + ^ l2+ ^ /2K2) \ t > 0, 

is a martingale with respect to <r(£ s : s < t). Define Q 1 ^ to be the prob- 
ability measure which has martingale density T K (t) with respect to 
on <t(£ s : s < t). Under Q^, £ is now a Brownian motion conditioned 
to stay in (0,K); cf. Knight [Kni69]. 

Using ideas in [HH09], we can use T K to construct a martingale with 
respect to T t = cr(X s ,s < t), the filtration generated by the P K - 
branching diffusion up to time t. For each u G Nt, write T^(t) = 
sm(irx u (t)/K)e^V + ^ 2 / 2+w2 / 2K ^\ t > 0. Define the process Z K = 
(Z K (t),t > 0) given by 

Z K (t) = e- (m - m T u (t) = e^ {t) - x{K)t sm(7rx u (t)/K), t > 0, 

u&N t u&N t 

where we set \{K) := (m — — fi 2 /2 — tt 2 /2K 2 . Then Z is a non- 
negative (P^, J-t)-martingale. For x G (0, K), we define a martingale 
change of measure on the probability space of the P^-branching diffu- 
sion via 

(2.2) 



z K {oy 



This change of measure induces the following spine construction for the 
path of X under ■ From the initial position x, we run a Q^f-diffusion, 
that is a Brownian motion conditioned to stay in (0, K), and we call 
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it a spine. At times of a Poisson process with rate mj3 we immigrate 
A independent copies of (X, P K ) rooted at the spatial position of the 
spine at this time. The number of immigrants A has the size-biased 
offspring distribution 

1 + k 

Qk = qk+i, k>0. 

m 

Let us remark that the change of measure with the martingale T K 
in (2.1) is equivalent to a Doob's /i-transform on L using h(x) = 
sm(7rx/K)e^ x . The infinitesimal generator L* of a Brownian motion 
conditioned to stay in (0, K) is therefore 

L* - 1 ^ I n / K ° d on (OK) 
K 2dx 2 tan.(TTx/Ko) dx 1 

with domain C 2 (0, K). We then call the Brownian motion conditioned 
to stay in (0, K) a L^-diffusion and note that it is positive recurrent 
with invariant density sm 2 (irx/K), for x £ (0, K). 



Remark 1. The martingale construction above applies more gen- 
erally to branching diffusions with spatially dependent branching mech- 
anism. Let us give a brief sketch as we will make use of it in Section 
4-2 (see again [HH09] for a full account). Suppose we have a branching 
diffusion Y on [0,K] with (spatially dependent) branching mechanism 
F(s,x), s £ [0,1], x £ (0,K) and set F'(x,l) := f s F(x, s)\ s=1 . Let 
T = (T(£),£ > 0) be a unit-mean martingale for the single particle 
motion and accordingly, for u £ N t , define T u (t) as the same object 
but with the associated single particle replaced by the particle position 
y u (t). Then 

(2.3) Z(t) = Vexpj- f F'{y u {s),l) ds)t u (t) 

defines a martingale with respect to cr(Y t , t > 0). Changing measure with 
the martingale Z induces a spine decomposition in which the spine has 
martingale density T(£) with respect to the law of the single particle 
motion in Y . 
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Let us continue with the study of the martingale Z . Since we as- 
sumed E(A\og + A) < oo, the following Proposition gives a necessary 
and sufficient condition for the L 1 (P^)-convergence of Z K . 

Proposition 1. Recall that \(K) = (m - l)/3 - /j, 2 /2 - n 2 /2K 2 
and let < x < K . 

(i) If X(K) > then the martingale Z K is L l (P^)- convergent and in 
particular uniformly integrable. 

(ii) IfX(K) < then ]im t ^ 00 Z K (t) = Pf-a.s. 

We refrain from giving the proof of Proposition 1 since it is a straight- 
forward adaptation of the proof of Theorem 13 in Kyprianou [Kyp04] 
which presents the /^-convergence result in the case of a branching 
Brownian motion with absorption at a space-time barrier, see also proof 
of Theorem 1 therein, as well as the proof in [Lyo97] and the proof of 
Theorem A in [LPP95]. 

We will now show that the martingale limit Z K (oo) is zero if and only 
if the process becomes extinct. 

Proposition 2. For x e (0,K), the events {Z K (oo) = 0} and 
{( < oo} agree Pfr-a.s. 

An essential idea in the proof of Proposition 2 is to embed the killed 
branching diffusion in a branching diffusion with killing on a larger 
strip. Let us introduce this procedure and some notation now as it will 
be used again in Section 4.2.2. 

Denote by Pi a ' b ^ the law under which X is our usual branching Brow- 
nian motion but with killing upon exiting the interval (a, b), where 
— oo<a<6<oo (and we simply write P% instead of Pj '^ in accor- 
dance with our previous notation). We denote by r u and o~ u the birth 
respectively death time of a particle u and write v < u if v is an ances- 
tor of u (u is considered to be an ancestor of itself), in accordance with 
the classical Ulam-Harris notation (see for instance [HH09], p. 290). For 
an e > 0, we choose a and b such that a < < e < b. Under Px a ' b \ we 
define 

^*l(o,e) = {ue N t : x v (s) G (0, e) Wv < u, r v <s< a v A t} , 

which is the set of particles u G Nt whose ancestors (not forgetting 
u itself) have not exited (0, e) up to time t. Now we can define the 
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restriction of X to (0, e) under Px by 

Xt\(o,e) = ^2 t - °- 

Then we conclude immediately that, for an initial position in (0, e), the 
restricted process X\^ 0e ^ = (x t \( 0f ^,t > OJ under Px a ' b ^ has the same 
law as (X t P£). 

Proof of Proposition 2. Clearly {( < oo} c {Z K (oo) = 0} 
and it remains to show that survival implies that Z(oo) is strictly 
positive. We consider the cases X(K) < and X(K) > separately. 
Assume X(K) < 0. Suppose for a contradiction that Z K (oo) = on 
survival. This requires the terms e^ Xu ^ sm(7rx u (t) / K) to vanish which 
can only happen if all particles move towards the killing boundaries 
and K. That is to say, for any e > 0, all particles leave the interval 
(e, K — e) eventually, and thus we may assume without loss of generality 
that the process survives in a small strip (0, e), for any e > 0. We will 
now lead this argument to a contradiction by showing that, for e small 
enough, the P^-branching diffusion, x G (0,e), will become extinct a.s. 
We embed the P e -branching diffusion in a P < -~ <5,e+<5 - ) -branching diffusion 
according to the previously described procedure. Now we choose 5 and e 
small enough such that A(e + 2<f) := (m- l)/3- y u 2 /2-7r 2 /2(e + 25) 2 < 0. 
Then, under p(- <5 > £ + ,5 ) ) the process 

Z (-*>'+f)(t) 

:= {e^M+< 5 )- A ( e+2 < 5 )' sin(7r(x u (t) + 5)/(e + 25))} , t > 0, 

ueN t 

is a martingale of the form in Proposition 1. Considering now the con- 
tribution coming from the particles in the set iVt|(o,e) on ly ; we nrs t 
note that survival of the P £ -branching diffusion ensures that this set is 
non-empty for any time t. Further, for particles u G Nt\(n e ), the terms 
e fi(x u (t)+S) sin^^x^t) -\- 8)/{e + 25)) are uniformly bounded from below 

by a constant c > and hence, under Px S ' e+S \ we get 

Z ( ^ +s \t)>cN t \ m e- x(e+25)t . 

Since we have chosen 5 and e such that A(e + 25) < 0, we now con- 
clude that Z(~ 5,e+5 \oc) = oo, Px~ S ' e+S ^-a..s. This is a contradiction since 
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is a positive martingale and therefore has a finite limit. Hence, 
for X(K) < 0, the martingale limit Z K (oo) cannot be zero on survival. 
Assume now \(K) > 0. Suppose for a contradiction that {( = oo} D 
{Z K (oo) = 0} is non empty and work on this event from now on. Now 
let z K {x) = P*{Z K (oc) = 0), for x £ (0, K). Define := 1 {z k {oo)=0} 
and set 

M t := E? (Mool^i) = J] z K (x u (t)), 

uGN t 

where the second equality follows from the branching Markov property. 
Then (M t ,t > 0) is a uniformly integrable P^-martingale with limit 
Mqo. Hence its limit on the event {( = oo}n{Z^(oo) = 0} is 1, P^f-a.s. 
This requires in turn that all particles x u (t),u £ N t move towards 
and K as t — > oo, since we know from Proposition 1 (i) that zk{x) < 1 
for x within (0, K). The previous part of this proof already showed that 
this leads to a contradiction. Thus, for X(K) > 0, the martingale limit 
cannot be zero on survival. This completes the proof. □ 

Proof of Theorem 1. Note that \(K) > if and only if [i < 
y/2(m-l)/3 and K > K . The result follows now immediately from 
Proposition 1 and Propositon 2. □ 

Remark 2. We can apply the same argument given for the process 
M in the proof of Proposition 2 to show that 

E«(l {C K <oo} \T t )= H (l-p K (x u (t))), t>0 

ueN t 

is a uniformly integrable product martingale which, followed by a clas- 
sical Feynman-Kac argument (cf. Champneys et al. [CHT^ 95]), gives 
that 1 — Pk(x) solves 

Lu + F(u) = on (0, K) 
(2.4) u(0) = u{K) = 1. 

We will show later in Remark 3 at the end of Section 3 that, if there 
exists a non trivial solution to (2.4), then it is unique. This will then 
again imply that {Z^oo) = 0} and < oo} agree P^-a.s. 
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3. Backbone decomposition via martingale changes of mea- 
sure. In this section we decompose the P^-branching diffusion into 
branching diffusions corresponding to the blue and red trees described 
in our intuitive picture in Section 1.1. The blue tree which consists of 
all genealogical lines of descent that will never become extinct will be 
shown to correspond to a branching diffusion which we will henceforth 
refer to as the backbone. Secondly, the red trees which contain all re- 
maining lines of descent will be shown to correspond to copies of the 
P K -branching diffusion conditioned on becoming extinct. 
The law of the branching diffusion corresponding to the coloured tree, 
P K , is defined by the law of X under P K and a subsequent deter- 
ministic colouring of the particles as described previously. Its natural 
filtration is T t '■= o~{T t , c(u) u£Nt } where c{u) is the colour of a particle 
u G N t . We say a particle u is blue if it has an infinite genealogical line 
of descent and we write c(u) = b, otherwise we say it is red and write 
c{u) = r. We have, for all t > 0, 



dP 



K 



n 



-{c{u)=b} + l{c(»=r} 



1 



and thus 
dP K 



dP? 



(3.1) 



E x( II ( 1 W«)=H + 1 {c(u)=r}") ?t 

ueN t 

e n p * r ( c(u) = cu V) 

j2 n n (}-PK{xu{t)) 

(cu)ueN t u€N t ,c u =b u£N t ,c u =r 



where (c u ) u£ N t is the set of all possible colourings of N t . In particular, 
for A G T t , we get 



P«(A;c(u) = c u \/uEN t \T t ] 



u£N t ,c u =b 



u£N t ,c u =r 



We can now easily derive the change of measure for the branching 
diffusion corresponding to the red tree. Let A G Tt, and write c(0) = r 
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for the event that the initial particle is red and thus 
P^):=P^|c(0) = r) = Pf(4««) = rW6tf t ) 



(3.2) 



Pf (c(0) = r) 

(iAn^ t (i-p*p^))) 



Clearly, conditioning the genealogical line of the initial particle to be- 
come extinct agrees with conditioning the whole process to become 
extinct and therefore the law of X under ~P R,K agrees with the law of 
X conditioned on extinction. The following Proposition characterises 
the process under ~p R,K . 

Proposition 3. For v e M a (0,K), define ~P R ' K via (3.2). Then 
(X, P R,K ) is a branching process with single particle motion charac- 
terised by the infinitesimal generator 

L « = l&-(^T^)i » <"■*>• 

for u G C 2 (0, K) with u(0+) = u(K-) = 0, and the branching activity 
is governed by the space- dependent branching mechanism 

Fg(s,y) = - 1 (F(s(l-p K (y))) - sF(l - p K (y))), 

1 -PkW 

for s G [0, 1] and y G (0, K) . 

Proof. The change of measure in (3.2) preserves the branching 
property in the following sense. Let v = Y^l=i be an initial configu- 
ration in (0, K) and A G Tt- Then 



EGU(i-p*0*)) 

fj^ Xi \ A 1 - Pk{xi) / 

(®klp2*)w 
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The process (X, P R > Rs ) is therefore completely characterised by its evo- 
lution up to the first branching time T. Let £ = {£ t , < t < T} denote 
the path of the initial particle up to time T and let if be a positive 
bounded measurable functional of this path. We begin with considering 
the case t < T. We have 



E£ K (H(t.,8<t);T>t) = E^H^sKt) ] PK f ) ;T>t) 

V l-p K (x) J 

(3.3) = e-^ K [H{^s<t)e' io ^k(M ds \ 

where P^-^ is defined by the change of measure 

and (Qt, t > 0) denotes the natural filtration of (£, P^). Thus the initial 
particle performs a P^'^-motion which is governed by the infinitesi- 
mal generator as given in the statement of the proposition. Taking 
H = 1 above, we see immediately that under ~P R,K the branching rate 
changes to 



(3.5) 



Pk{v) 



for y G (0, K). 

It remains to identify the offspring distribution and we therefore study 
the process at its first branching time T. We get 

E£*(ff(&, s < T); N T = k- T E dt) 

- ^{ {1 ;^ T m„s<nT,^N T ^ 



E 
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We see that, in addition to the change in the motion and the branching 
rate, the offspring distribution under P^-^ becomes {q R , k > 0} where 

(3.6) q R (y) = ^ R {y)y l q k {l - VK {y)) k ~ l , k > 0. 

A simple computation shows that F^(s, y) = f3 R (y)(J2k>o Qk(v) sk ~ s ) 
takes the desired form. □ 

Similarly to the above reasoning, we obtain the law of the dressed 
backbone, that is a backbone with immigration of P R,/< -branching dif- 
fusions, by conditioning on the first particle being blue. Thus 

P°' K (A) := Pf(A|c(0) = 6) 

P^(A; c(u) = b for at least one u £ N t ) 



(3.7) 



Pf (c(0) = b) 
p K (x) 



Then (X, P D > K ) certainly agrees with (X, P ) conditioned on survival. 
Let us characterise the evolution under P® ,K ■ We use the previous 
notation and in addition let r = T A VQ t K) denote the death time of the 
initial particle, where T7o,in is the first time this particle exits (0,K). 
Then 



E 



D,K / 



(H(Z s ,s<t)-T>t) = e-^(H(C s ,s<t)^l,T i0)K) >t) 

\ Pk{x) / 

(3.8) = e~^ K [H (&, s < t) ds ) , 

where V B,K is defined by the change of measure, for t > 0, 



L { r (0,K)>*}' 



(3.9) 

and > 0) is again the natural filtration of (£, Pf"). Thus, setting 

A*) = 



N F(l -p^(x)) - /3p^(x) 



(3.10) = /j 1 -^ 1 -^^ , forx£(0,i,), 

Pk{x) 
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we see that (3.8) simplifies to 

Ef s < t),r > t) = EZ> k (h(£ s ,s< t)e-fi pD M ds 
(3.11) 

We deduce from this that, under P D,K , the motion of the initial parti- 
cle is given by the change of measure in (3.9) and it branches at rate 
(3 D {-) as in (3.10). 

It remains to specify the offspring distribution. We begin with the ex- 
pression in (3.7) and then use (3.9) and (3.10) to get 

Ef s<T);T edt;N T = k) 

\k 



V vk lx) / 



p K {X 

v P d \Kt) Pk[£t) 

Again this reveals the evolution of the initial particle as described above 
and we further see that the offspring distribution of the initial particle 
under P D ' K is given by {q^, k > 0} where 

up to the normalising constant (3(f3 D (x))~ 1 . We note that qo(x) = 
for all x G (0, K) which we expected to see since (X, P D < K ) is equal 
in law to (X, P K ) conditioned on survival. However, we have so far 
neglected the fact that the initial particle can give birth to particles 
of the same type, i.e. blue particles (referred to as branching), and red 
particles which evolve as under P R ' K (referred to as immigration). We 
will split up the rate (3 D and the offspring distribution into terms 
corresponding to branching respectively immigration. Firstly, note that 
we can decompose the rate (3 D into 

k>2 n>k ^ ' 

+ pY J q n n{i-VK{y)) n - 1 

n>l 

(3.12) =: (3 B (y) + (3 I (y). 
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Then (3 1 is the rate at which the initial particle gives birth to one blue 
particle and a random number of (red) immigrants (immigration rate) 
while j3 B is the rate at which the initial particle gives birth to at least 
two particles of the blue type and a random number of (red) immigrants 
occur (branching rate of the branching diffusion corresponding to the 
blue tree). We can now rewrite the offspring distribution as 

1 - {1 - Pk ( X ))*t 

Pk{x) 

(3.13) = q k Y,( k )pK(x) i - 1 (l-p K (x)) k - i 

i=2 ^ ' 

(3.14) +q k (l-p K (x)) k -\ k>\. 

Then the term in (3.13) gives, up to normalisation, the probability that 
the initial particle branches into i particles of its type and, at the same 
branching time, k — i particles immigrate. The term in (3.14) is the 
the probability that k — 1 immigrants occur, again up to a normalising 
constant. 

Note that (X, P D > K ) inherits the branching Markov property from 
(X,P K ) by (3.7) in a similar spirit to the case of (X,P R ' K ) (cf. the 
proof of Proposition 3). Thus the description of the initial particle also 
characterises the evolution of all particles of the blue type and together 
with the characterisation of the immigrating P R '^-branching diffusions 
in Proposition 3 we have completely characterised the evolution of X 
under P D ' K . The following result is now an immediate consequence. 

Theorem 5 (The dressed Backbone). Let K > K and x 6 (0, K) . 
The process (X, P^ ,K ) evolves as follows. 

(i) From x, we run a F^' K -diffusion which dies at rate j3 B , 

(ii) at the space-time position of its death, it is replaced by A B parti- 
cles where A B is distributed according to the probabilities, 

Qk(y) = P ^(yr^^QpKiyf^ii-PKiy))^, 

n>k ^ ' 

(3.15) 

fork>2 and y 6 (0, K) . 
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(Hi) Each of the offspring particles repeats its parent's stochastic be- 
haviour. 

(iv) Conditionally on the branching diffusion, say X B , generated by 
steps (i) - (Hi), we have the following. 

• ( Continuous immigration) Along the trajectories of each par- 
ticle in X B , an immigration with n > 1 immigrants occurs 
at rate 

#(y) = 0q n+1 (n + l)(l-p K (y)) n , y E (0,K). 
(3.16) 

• (Branch point immigration) At a branch point of X s with 
k>2 particles, we see an immigration of n > immigrants 
with probability 

qi(y) = q n+k ( U + % K {yf-\l - p K (y)) n , y E (0, K). 



k 



(3.17) 



Each immigrant initiates an independent copy of (X, ~P R,K ) from 
the space-time position of its birth. 

The following Theorem is a more precise version of Theorem 2, stated 
in Section 1.1. 

THEOREM 6 (Backbone decomposition). Let K > K and v E 
M a (0,K) such that v = £)" =1 5 X . with x { E (0,K), n > 1. For t > 0, 
we can define 



dP 



K 



dP? 

(3.18 



n k rlT> D > K 

= E E Up^ipr nd-pfe)) 

Tt k=0 (x!,...,x k ) i=l Xi 



dP 



R.K 



dP* 

Tt j=k+l X J 



Ti 



where the second sum above is taken over all k-tuples ofx\, ...,x n . How- 
ever note that the right-hand side of (3.18) is equal to 1 and thus, triv- 
ially, on the filtration (J"t)t>o, {X, P^) is Markovian and equal in law 
to (X,P?). 
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Proof. The change of measure is just a restatement of (3.1) and 
the result follows from Proposition 3 and Theorem 5. □ 

Intuitively speaking, we can describe the evolution under P® ,K and 
thus also under Pff as follows. Independently for each initial particle 
Xi, we flip a coin with probability p(xi) of 'heads'. If it lands 'heads', 
we initiate a copy of (X, P^ K ) and otherwise we initiate a copy of 

(*,p5*)- 

Corollary 1. Given the number of particles of (X, P^) and their 
positions, say x±, x n for some n G N ; at a fixed time t, the number of 
particles of Xf is the number of successes in a sequence of n indepen- 
dent Bernoulli trials each with success probability pk{xi), ■■■,PK{x n )- 

We refer to the branching diffusion generated by steps (i)-(iii) of 
Theorem 6 above as the backbone. Its law can be characterised as 
follows. 

Proposition 4 (The backbone). For v e M a (0, K) such that v = 
Y^i=i^xi with Xi G (0, K), n > 1, we define the measure P B,K via the 
following change of measure. For t > 0, 



dP u 



B,K 



dP, K 



-j-r p K {Xy{(Tv At)) 

11 PK (x v (r v )) {4< W 



{/•(TvAt \ 
J F'{l-p K {x v {s))) + (3 dsj 



_qAjM°vm B {x v {a v )))^ 

where % is the set of all particles v G T with r v < t and v is in Tt- if, in 
addition, a v < t . As usual, r v and a v are the birth respectively death 
times, t? q K ^ is the first exit time from (0, K) and A v is the random 
number of offspring of a particle v G %- ■ 
The branching diffusion (X, P„ ,K ) has infinitesimal generator 
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defined for all u G C 2 (0, K), and space- dependent branching mechanism 

Fk(s, V) = — U (F(sp K (y) + (1 - PK (y))) - (1 - s)F(l - p K (y))) , 
PK{y) 

for s G [0, 1] and y G (0, K). The process (X, P B ' ) evolves according 
to the steps (i)-(iii) of Theorem 6. 

Proof. First note that the motion under F B ' K , given by the change 
of measure in (3.9), is governed by the infinitesimal generator as 
given in the statement. A simple computation also shows that F^(s, y) = 
P B (y)(J2k>2<lk( x ) sk ~ s ) with (5 B and q B as in (3.12) and (3.15) gives 
the desired form. The result then follows from rewriting the change of 
measure up to the first branching time T as 



dP 



B,K 



dP? 



Pk(x) [ J p K (i 
xi/3 B (£T))exp j- £ (3 B ^ S ) - (3 ds 



ex P \ ~ I — t^jt\ — ds r 1 {*<-(o,x ) } 

X 



noting that the first line on the right-hand side accounts for the change 
of motion, the first term in the second line for the change in the branch- 
ing rate and the last term in the second line for the change in the 
offspring distribution. □ 

Remark 3. As promised earlier, with the help of Corollary 1, we 
can show that, if (2.4) has a non-trivial solution, then it is unique. 
Assume gxix) is a non-trivial solution to (2.4). By a Feynman-Kac 
argument (cf. Champneys et al. [CHT^ 95]), it follows that 

M K (t) = J] g K (x u (t)), t > 0, 
ueN t 

is a -product martingale. Since M K is uniformly integrable, its limit 
M K (oc) exits P^-a.s. On the event of extinction, M K (oo) = 1. On the 
event of survival, we have 

(3.19) M K (t) = J] g K (x u (t)) < J] g K (x B (t)), 
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where N B is the set and x u (t) are the spatial positions of the par- 
ticles in Xf . Clearly \N t \ B — > oo as t —> oo since each particle in 
X B is replaced by at least two offspring and there is no killing. Fur- 
ther particles in X B perform an ergodic motion and it is therefore 
not possible that liminf g(x u (t)) tends to 1. Thus the right-hand side 
of (3.19) tends to and we conclude that M K (oo) = l^ <00 y Hence 
gx^x) = E% (M K (oo)) = P^(( < °°) which implies uniqueness. 
In particular we may conclude that (2.4) has a non-trivial solution if 
and only if ' /i < \fl(m —1)13 and K > Kq . 

4. Proof of Theorem 3. We break up Theorem 3 into two parts 
which will be proved in the subsequent sections. 

Proposition 5. Uniformly for all x e (0, K ), 

Pk(x) ~ ck sm(iix/K )e tlx , as K I K , 

where ck is independent of x and ck I as K I K . 



Proposition 6. The constant ck in Proposition 5 satisfies 
(4.1) c K ~ (^-^) 12(m _ 1)/37r ^ 3(e ^ + 1) asKlK . 

Theorem 3 then follows by defining Ck to be the expression on the 
left-hand side in (4.1). 

We will provide entirely probabilistic proofs of the results above. We 
remark that, although it would take some effort to make rigorous, it is 
also possible to recover the asymptotics of px and the explicit constant 
Ck in an analytic approach using a careful asymptotic expansion of 
the non-linear ODE Lu + F(u) = with u(0) = u(K) = 1, as shown 
to us by B. Derrida. 



4.1. Proof of Proposition 5. We begin with a preliminary result 
which ensures that the survival probability px is right-continuous at 

Lemma 1. Let x e (0,Ko). Then lim K i Ko p K (x) = 0. 
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Proof. We fix x £ (0, Kq) throughout the proof and consider Pk{x) 
as a function in K. For a fixed t > 0, let us define the probabil- 
ity px{x,t) := P^survival in (0,K) up to time t). By monotonicity 
of measures we have \im K i Ko p K (x, t) = px (x,t). Now we can write 
Pk{x) = ini t>0 pK(x,t). Hence Pk(x) is the infimum of a sequence of 
functions which are continuous at Kq and thus upper semicontinuous 
at K , that is 

lim sup ^(x) < Pk (x). 

KIK 

Furthermore, Pk{x) is decreasing as K 1 Kq and bounded, so the right 
limit exists and 



Pk (x) < lim p K (x). 
A4.A0 

Combining the two inequalities above we obtain right-continuity of 
Pk{x) at Kq. By Theorem 1, Pk ( x ) — and so we have lim/qj^ Pk(x) = 
0. 

□ 



The following lemma is the essential part in the proof of Proposition 
5. Its proof is guided by the ideas in Aidekon and Harris [AH]. 

Lemma 2. Let y £ (0, Kq). Then we have 
(4.2) lim ^ = S ^ x/K °]e^-y\ 

uniformly for all x £ (0, Kq). 

Proof. Fix y £ (0,Kq). We begin with showing that the asymp- 
totics hold uniformly for all x £ (0,y). 

Let T denote the set of labels of particles realised in (X,P^). Define 
T y as the set containing all particles which are the first ones in their 
genealogical line to exit the strip (0,y), i.e. 

T y = {u £ T : 3s £ [t u , a u ] s.t. x u (s) <£ (0, y) 

and x v {r) £ (0,y) for all v < u,r £ [t„, <7„]}, 
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where v < u means that v is a strict ancestor of u. Further, for u G T y 
denote by T y the first exit time of u from (0, y). The random set T y is 
a stopping line in the sense of Biggins and Kyprianou [BK97] (see also 
Chauvin [Cha91] which uses a slightly different terminology though). 
Since y G (0, K ) the width of the strip (0, y) is subcritical and hence, 
for any initial position x G {0,y), all particles will exit it eventually 
which ensures that T y is a dissecting stopping line. Now let \T y \ be the 
number of particles which are the first ones in their line of descent to 
hit y, which can be written as 

\ T y\= ^2 1 {MT y u )=y}- 

ueTy 

Recall from Remark 2 that QXteiV ~ PK(x u {t))),t > 0) is a P^- 
martingale. Since T y is dissecting it follows from [Cha91] that we can 
stop the martingale at T y and obtain, for x G (0, y), 



l-p K (x) = E« [ J] l-p K (x u (T-)) ) = EX((l-p K (y))\ T y\) 
(4.3) 



k «6T„ 



where we have used that the process started at zero becomes extinct 
immediately, i.e. Pk{0) — 0. Further \T y \ has the same distribution 
under and P^° since we consider particles stopped at level y below 
K Q and thus we can replace by E^° on the right-hand side above. 
Now, using first (4.3) and then the geometric sum Y^j=o aJ = T^T> we 
get 

^ = E ''° { l ~i ( \7 K{ T^ ) = E -° I ^ < J - ^ 
pkw \ i-{i-p K {y)) J 

(4.4) 

The sum on the right-hand side is dominated by \T y \ which does not 
depend on K and has finite expectation (which will shortly be shown 
below). We can therefore apply the Dominated convergence theorem 
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to the right-hand side in (4.4) and we conclude that 

(\Ty\~l 



3=0 



\Ty\-l 



(4-5) = E«° I g Km (1 -pK(2/)) J j = ^°(|T W |), 

where the convergence holds point-wise in x G (0, y). In order to 
get uniform convergence we observe the following. We set ip(x, K) = 

E^° (^2^=0 X (l — PK:(y)y\ for x G [0, y] (with the convention that 

the P^-branching diffusion becomes extinct immediately for initial po- 
sition x = respectively stopped for x = y) and denote by <f(x) = 
.E^ (|T y |) its point-wise limit. Since l—pxiy) < l—PK'(y), for K > K', 
we have ip(x,K) < (p(x,K') and thus, for any x G [0, y], the sequence 
<p(x, K) is monotone increasing as K 1 Kq. Moreover the functions 
<p(x, K) and (p(x) are continuous in x, for any K. In conclusion, we 
have an increasing sequence of continuous functions on a compact set 
with a continuous point-wise limit and therefore the convergence also 
holds uniformly in x G [0,y] (see e.g. [Rud76], Theorem 7.13). 
Combining (4.4) and (4.5) and the uniformity argument, we arrive at 

(4.6) Km = E«>{\T V \), 

KiK p K (y) 

where, for fixed y, the convergence holds uniformly in x G (0, y). Now 
let Tg := inf-jr > : £ t G (0, y)} be the first time a Brownian motion £ 
with drift —y, exists the interval (0, y). Since T y is dissecting it follows 
from Theorem 6 in [Kyp04] that we can apply the Many-to-one Lemma 
(see e.g. [HH09] Theorem 8.5) for the stopping line T y . This gives 




sm(nx/K )e^ {m -i)^ v , 

sin(vr^/K )e^ + ^/ 2+ - 2 / 2 ^)^ ' ^ =V) 

sm(iry/K ) iv 
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where we have used that (m — l)/3 — \i jl — tt 2 /2Kq = (and (Q£° is 
used as an expectation operator). Under Q^°, £ will never hit since 
it is conditioned to stay in (0,K ). However as £ is positive recurrent 
it will eventually cross y and therefore Q^°(^ y = y) = 1- This proves 
our earlier claim that \T y \ has finite expectation and together with (4.6) 
it completes the argument. 

For uniformity for all x G (0,K), it remains to show that (4.2) also 
holds uniformly for x G (y,K ). Instead of approaching criticality by 
taking the limit in K we can now fix a K > K and consider a (su- 
percritical) strip (z, K) and let z t z where z := K — K . Denote by 
P( z k)( x +z) the probability of survival in the strip (z, K) when starting 
from x + z. We then have 

lim^M = lim^ )(X + " ) . 

KIKo p K (y) ztzo P( z ,K)(y + Z) 

Hence (4.2) is equivalent to showing that, uniformly for x G (y,K Q ), 

*tzo P( Z , K ) (y + z) sm(ny/ K ) 

Then consider the stopping line containing all particles which exit the 
strip (y + z, K) and accordingly the set of particles which are the first in 
their genealogical line to exit (y + z, K) at y. Noting that the latter has 
the same law under P%+ z and P^+ z , we can then repeat the argument 
in the first part. □ 

Proof of Proposition 5. Choose aye (0, K ). Then an appli- 
cation of Lemma 2 gives, as K j. K , 

, \ , \Vk(x) , x sin(7Tx/i^ ) u(x-v) I its\ ux 

Pk(x) = p K (y) — K ~ PK(y)— ( — T^ e = c K sm(7rx K )e^ x , 
p K (y) Bm(iry/K ) 

uniformly for all x G (0, K ), where ck '■= ih^^^T e ~ ra - By Proposi- 
tion 1, ck 4- as K 4, K which completes the proof. □ 



4.2. Proof of Proposition 6. In this section we will present the proof 
of Proposition 6 which gives an explicit expression for the constant ck 
appearing in the asymptotics for the survival probability in Proposi- 
tion 5. We outline a heuristic derivation of the explicit constant ck 
in Proposition 6 which will give the intuition for the rigorous proofs 
presented subsequently. 
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4.2.1. Heuristic argument. We break up the heuristic argument for 
Proposition 6 into four steps which we will refer back to in the rigorous 
proof in the subsequent section. 

Step (i) (The growth rate of the backbone) Consider a process Y B = 
{Yt B >t — 0) performing the single particle motion of the backbone, 
that is according to the infinitesimal generator which is given in 
Proposition 4 as 

with domain C 2 (0, K). Let IT^ be the invariant density for L^, i.e. the 
positive solution of L^U 1 ^- = where is the formal adjoint of L B . 
Then 

nf (x) oc p K {x) 2 e- 2 >* x , x e (0, K). 

For t > , we define T(t,A) = f Q l{Y 3 B eA\d> s i A C [0, K], to be the 
occupation time up to t of Y B in A. Then large deviation theory 
suggests that the probability that the measure * -1 r(t, •) is 'close' to 
J K l { . } (y)f(y)Il^(y)dy should be roughly 



e XV {t£ L B f(y) f{y)U*{y) dy 



Now, as each particle in the backbone moves according to L^, we guess 
that the expected number of particles at time t with occupation density 
like / 2 n^ is very roughly 

(4.7) exp jt J (L| + F#'(l, y))f(y) f(y)U B (y) dy 
where 

x) := ^-F B (s, x)U = (m - + g^Z**^ , x G (0 , X). 
as Pk{x) 

The expected growth rate of the blue tree is^iven by maximising the 
integral appearing in (4.7) over all / with j Q f 2 (x)U B : (x)dx = 1. We 
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can compute this optimal function /* explicitly as the normalised eigen- 
function corresponding to the largest eigenvalue A where 

(4.8) (L B K + F*'(l,x))r(x) = \r(x) m(0,K), 

and we find that, in fact, A = X(K) = (m — l)/3 — /i 2 /2 — ix 2 /2K 2 and 

(4.9) rtooc ^/fV , xe(0,K), 

up to a normalising constant. Then we find the 'optimal' occupation 
density as 

<'*(x) := (/* (x)) 2 nf (x) = A sin 2 (vrx/K), x G (0, K). 

In summary, we guess that the expected growth rate of the number of 
particles in the blue tree is X(K) and that 

(4-10) X(K)= f K (L*+F*\l,y))r(y) r(y)U B K (y) dy. 
Jo 

We would anticipate that the a.s. growth rate is also X(K) in agreement 
with the expected growth rate. 

Step (ii)(Upper bound on X(K)) The term L B : f*f* is non-positive as 
it represents the cost of spending time like (/*) 2 IT^-, hence omitting it 
will give an upper bound for X(K), that is 

-K 



f F B '(l iy )U B K (y)dy < X(K) 
Jo 



Step (iii)(Lower bound on X(K)) By taking / = 1 in (4.10), we get a 
lower bound on X(K) since /* maximizes the expression in (4.7). Thus 



X{K) < f K F*'(l,y)IL%'(y)dy. 
Jo 

Step (iv) (Asymptotics) By Theorem 3, Pk( x ) ~ c k sm(7rx/K )e tJ,x , as 
K l K , and we can easily deduce that Hk(x) ~ Il^' o *(x), as K I K . 
We will make rigorous later that F B '(1, x) ~ (m — l)/3c^ sin(7ra;/ 'K )e' J " x 
as K I K . Our conjecture is therefore that 

X(K) ~ c K -f / sm 3 (iry/K )eM dy, as K | K . 
A o Jo 
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Since we can calculate the integral explicitly this gives an exact asymp- 
totic for ck which agrees with the one given in Proposition 6 and The- 
orem 3. 

4.2.2. Proof of Proposition 6. Let us now come to the rigorous 
proof. Recall that the backbone (X, P B > K ) is the process constructed in 
steps (i)-(iii) in Theorem 6, which was further characterised in Propo- 
sition 4. First, we need to confirm the conjecture that the number of 
particles in (X,P^ ,K ) grows at rate X(K). 

Proposition 7. For x e (0, K), 

lim -log \N t \ = X(K), P^ K -a.s. 

t— too t 

Step (i) of the heuristic suggests that the growth rate of the backbone 
is X(K) and, moreover, that it can be expressed as (4.10). The idea of 
the rigorous proof of Proposition 7 is now to construct a martingale 
of the form in (2.3) with T built from /* in (4.9). We then find upper 
and lower bounds for this martingale which will, in turn, give bounds 
on the growth of the number of particles. Before we do this we prove 
an auxiliary result on the boundedness of /*. 

Lemma 3. The function f* is uniformly bounded in (0,K). 

Proof. The function /* is continuous in (0, K) and it is there- 
fore sufficient to show that limsup x i f*(x) and lim sup^^ /* (x) are 
bounded. 

An application of L'Hopitals rule gives 



(4.11) lim .sin(W^K* 



3^(e- 2 ^ - 1) 

To conclude that lim sup f*(x) < oo, it therefore suffices to show 
that there exists a constant c > such that 

c (1 — e 2 ^ x ) < p K (x), for all x sufficiently close to zero. 

— 2/i 

By Remark 2, (Iluejv 0- ~ PK(x u (i))),t > 0) is a P^-martingale and 
it follows then by a standard Feynman-Kac argument that 1 — px{x) 
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satisfies 



l-p K (x) = 1 + E; 



X 



K 



f 

Jo 



T (0,K) 



F(l-p K (Q) ds, xe(0,K), 



where T(o,k) is the first time £ exists the interval (0, K). To compute the 
expectation above we use the potential density of £, see e.g. Theorem 
8.7 in [Kyp06], and we get 



Since F(s) < for < s < 1, the first integral in the last equality on the 
right-hand side of (4.12) is strictly negative. Regarding boundedness of 
this integral it is clear that the integrand is bounded for y near K. By 
an application of L'Hopitals rule it follows that the integrand is also 
bounded near 0. Hence we can set 



With the second integral in the last equality on the right-hand side of 
(4.12) being non-positive, for x close to 0, we get 

Pk(x) > 2c (e~ 2fix — 1) for all x sufficiently close to zero. 



which, by (4.11), gives the desired result. 

To establish boundedness as x approaches K, we observe that px(x) = 
Pk{K — x), where px denotes the survival probability for a branching 
diffusion which evolves as under but with positive drift fi. Similarly 
to the previous argument we can then show that there exists a constant 
c > such that cpx{K — x) > sin(7rx/A')e MX , for x sufficiently close to 
K, which finishes the proof. □ 

Proof of Proposition 7. We proof the upper bound by contra- 
diction. Recall the embedding procedure described in Section 2. Choose 



(4.12) 





dy > 0. 



-2/i 
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e small enough such that A(A + 2e) > 0. Choose a 5 > and suppose 
that there exists an increasing (random) sequence t n , n = 1, 2, which 
tends to infinity, such that log A 4n is bigger than (A(A + 2e) + S)t n , for 
any n E N under P K . Then, under p( e >- K '+ e ) ) 

Z {e ^ +e) (t) := sin(7r(x u (t) + e)/(KT + 2e ))e- A( ^ +2e V (x " w+e) 

> l^l(o,^) I e- x{K+2 ^ x sin(vre/(A' + 2e)). 

Since we assumed that, along the sequence t n , n = 1, 2, the number 
of particles \N tn \r K \ | is bounded from below by exp{(X(K+2e)+5)t n }, 
the right-hand side above tends to infinity along this sequence as n — > 
oo which contradicts Proposition 1. As we can take e and 5 arbitrary 
small we obtain limsup i _ il00 (A(A')t)~ 1 log \N t \ < 1, under P K . By the 
thinning argument in Corollary 1, we immediately get that A (A) is also 
an upper bound for the growth rate of \N t \ under P B,K . 
For the lower bound, as alluded to above, we begin with constructing 
a P B '^-martingale of the form (2.3). Since /* satisfies (Lg + P#'(l) - 
A (A))/* = 0, it follows by an application of Ito's formula that 

is a martingale with respect to cr(£ 4 ,t > 0), where (^,F B,K ) is an L^- 
diffusion. Appealing to the discussion in Remark 1 we then see that 

M r (t)= ^/*(* u (t))e-W t>0, 

uGN t 

is a Pf '^--martingale. 

The proof of L 1 (P^ ,/< )-convergence of Mj* follows by a classical spine 
decomposition argument (cf. the proof of Theorem 1 in [Kyp04], as 
well as the proofs in [Lyo97] and [LPP95]) and is therefore omitted. 
/^(Pf '^-convergence implies then that P B ' K (M f *(oo) > 0) > 0. Now 
set g(x) := ¥ B ' K {M f *{oo) = 0), for x E (0, A). Then the product 

* g (t) = II 9Mt)), t > 0, 

ueN t 

is a P^' x -martingale with almost sure limit l{M r (oo)=o} (cf. proof of 
Proposition 1). Therefore we have 

g(x) = V B ' K (^(t)) < Ef ^(gm, for all x E (0, A). 
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Hence we conclude that the process (g(£t), t > 0) is a [0, l]-valued ' - 
sub martingale and it converges P^'^-a.s. to a limit . However £ is 
positive recurrent under ¥ B,K and thus g{£,t) can only converge if it is 
constant, hence g(x) = for all x G (0, K). Since < g < 1 we then 
have g^ G [0, 1]. Assume now that g^ G [0, 1). Note that, under ,K , 
\N t \ tends to infinity as t — )■ 00 since each particle in (X, is re- 

placed by at least two offspring when it dies and there is no killing. Thus 
we get ir 9 (t) 0, Pf'-^-a.s. and therefore g(x) = E^' K (n 9 (00)) = 0. In 
conclusion, g is identical to either or 1. But we already know that the 
martingale limit Mj*(oo) is strictly positive with positive probability 
and consequently, g(x) = P^' K (M f ,(oo) = 0) = 0, for all x G {0,K). 
We may now conclude that lim inf log Mf . (t)/\(K)t > Pf'^-a.s. 
We look for an upper bound on Mf*(t) which will, in turn, provide a 
lower bound on \N t \ under P% ,K . By Lemma 3, /* is bounded by a 
constant c > in (0, K). Thus, under P^ ,K , for t > 0, 



M f *(t) <c\N t \ e 
and we see that, P^' K -a.s., 



-X(K)t 



l iminf ffl > ^ogMrV-logc + xmt 

t^oo X(K)t - t^oo X(K)t 
which completes the proof. □ 

In step (ii) of the heuristic we claimed that (L^f*)f* is non-positive 
to get an upper bound on X(K) which is essentially what we will now 
do. Recall that the invariant density for the infinitesimal generator Lf^- 
respectively the optimal occupation density of (X, P^' K ) are given by 



Jo Pk{ z ) 2 g~ 2 ^ z dz 
(4.13) and Xl B /{y) := (f* {y)f Hf (y) = A sm 2 (ny/K). 

Lemma 4. For K > K , we have 

-K 



X{K)< I F»\l,y)Il B K '*(y)dy. 
Jo 
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Proof. We have /* (f*(y) ) 2 ng(y) dy = J* 2/K sm 2 (7ix/K) dx = 
1. Then multiplying by \(K) gives 

K K ) = C KK)r(y)r(y)U B K (y) dy 
Jo 

Recall that /* is given by (4.9) and satisfies the ODE in (4.8). Thus we 
can replace the term X(K)f* above by (L B +F^ (1, y))f*(y). Therefore 

K K ) = ^^(r(«)r-(A*-^^)(/w)r(y)n?(y)dy 

(4-14) + f K F*'(l,y)(r(y)) 2 n B K (y)dy. 

Jo 

Noting that = (/*(y)) 2 nf-(y), the result then follows if we can 

show that the first integral in (4.14) is non-positive. 
We use integration by parts for the first term in the first integral in 
(4.14) to get 

K \iny))"rm B K{y) ^ 
= l(i(r(y)yrm B K (y)} K 

K (F(y))' (CT (l/))H£ (y) dy + r(y)(U B K (y))' dy) ) . 
jo 7 

(4.15) 

We want to show that the first term on the right-hand side above is 
zero. By Lemma 3, /* takes a finite value at and K and hence it 
suffices to show that (/*)TI^ evaluated at and K is zero. By simply 
differentiating /* and recalling that Ilf-(y) oc PK(y) 2 e~ 2fiy we get 



(r(y)yn B (y) oc e-^^sm(7ry/K) + -cos(7ry/K))p K (y) 

-sm(iry/K)p' K (y) 

Differentiating both sides of equation (4.12) with respect to x, it is eas- 
ily seen that p' K (x) is bounded for all x G [0, K\. Therefore (f*{y) )'H B : (y) 
is equal to at and K and thus the first term on the right-hand side 
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of (4.15) vanishes. 

The first integral in (4.14) now becomes 

£ QcrGOr - (a* - j^N^y))') r(v)n£G/) 

* (((/*(y))') 2 - (r (i/))T(y)(nf (y))') n£(„) dy 
(4.16) - A/, - (y)Cf (y))'n£ (y) 

Differentiating nf gives (ITf )' = -2(/i - ^)n£ . Thus, in the right- 
hand side of (4.16), the second term in the first integral cancels with 
the second integral and we arrive at 

£ Q(f (y))" -0*- /*(y)nf (y) d y 

= -\ J ((r(y))') 2 nf (y) dy, 

which is less than or equal to zero and the proof is complete. □ 

Step (iii) of the heuristic claims that we can lower bound X(K) by 
replacing /* in (4.9) with / = 1. This suggests to modify the martingale 
argument in the proof of Proposition 7 using a martingale of the form 
(2.3) with T built from the constant function 1. 

Lemma 5. For x e (0, K), 

\{K)> I <(i,y)nf(y) dy. 

Jo 

Proof. The constant process T = (T(£) = 1, t > 0) is a trivial mar- 
tingale with respect to er(£t, i > 0), where (£, ~P B)K ) is a /^-diffusion. 
Thus according to Remark 1, the process Mi = (M\(t),t > 0) defined 
by 

«eJV t B 
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is a Pf '^-martingale. L 1 (P^ ,i ^)-convergence and uniform integrability 
of Mi again follow by a spine decompisition argument in the manner 
of the proof of Theorem 1 in [Kyp04]. L 1 (P^'^)-convergence implies 
that P^ ,K (Mi(oo) > 0) > and repeating the argument in the proof 
of Proposition 7 we immediately see that Pf'^(Mi(oo) > 0) = 1. 
Therefore, we conclude that liminf^oo log M 1 (t) / \(K)t > 0, Pf'^-a.s. 
On the filtration Tt = (?(X S , s <t) we define a change of measure by 

dP^ K MAt) 

t > 0. 



M^oy 



Since the martingale Mi is of the form in (2.3), it induces as spine 
decomposition which, according to Remark 1, reads as follows. Under 
P^' K , the spine £ is an L^-diffusion and along its path we immigrate 
independent copies of the P^'^-branching diffusion (we do not need to 
specify the rate of immigration and the distribution of the number of 
immigrants since they will not be relevant). 
Next, fix an e > and define, for t > and each u G Nf, the set 



At 



r-K 



F» (1, x u (s)) ds - / Fg (1, y) IL%(y) dy 



< e 



and consider the process we obtain from Mi by considering the particles 
in Af only, that is 

Mi(t) = 1 Ar e-ti F K(^(s)) <fa t > Q 

Let A\ be the event we get if we simply replace x u {t) by the spine 
process £ t in the definition of Af. Since (£, p Ml '^) has invariant density 
II^- we have l A t — > 1 P^ 1 ' -a.s. By Theorem 1 in [HR], Mi therefore 
has the same limit as Mi under P^ 1,K , and moreover, since Mi is 
uniformly integrable, this also holds true under P^ ,K . In particular we 
have 

logM^) = logi^) 

t^oo A(K)t fr-*x> A(K)t _ 

For t > 0, we now get an upper bound for Mi(i) under Pf'^ by 
Mi(t) < j^ie-*^^'^) n f GO d v-\ 



36 HARRIS, S.C., HESSE, M. AND KYPRIANOU, A.E. 

Consequently, P^' -a.s., 

r • f log |iV t | 

Iim mi 



t— >oo 



> lim inf 



t^Fi'ihy^iy) dy 
log M^+t (/* F£' {l,y)Xl B K {y) - e 



> 



*(/*^Uv)nf(i/)di/ 

f?F*{l,y)IL%(y)dy-e 



f K F*'(l,y)n*(y) dy 
and taking e J, gives the result. □ 

Proof of Proposition 6. By Lemma 4 and 5, we get the follow- 
ing bounds on X(K) 

[ K Fi\l,y)Jl B K {y) dy < \{K) < f* F*'(l, y)IL%'(y) dy, 
Jo Jo 

(4.17) 

where IT^ and IT^'* were defined in (4.13). By Proposition 5, we have, 
as K I K , 

Jo Pk{z) e~ 2flz dz A o 

(4.18) 

where we have used that the asymptotics in Proposition 5 hold uni- 
formly to deal with the integral in the denominator of the second term 
in (4.18). The uniformity in Proposition 5 also ensures that (4.18) holds 
uniformly for all y G (0, K ). Further 

F( S ) /3(En> 2 gn^-l) ^ n _ x 

Iim — = lim = = limn > q n ns = [m — 1 p, 

v ' n>2 

where we applied L'Hopitals rule in the second equality above. Then, 
together with Proposition 5, as K I K , 

Fi'(l,y) = {m - 1)0 ^ -'*<»)) 

PkW 

(4.19) ~ (m-l)(3p K (y)^(m-l)(3c K sm(7ry/K )e^. 
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Note that f ^~pM = fi-^fa)) _ Convexity of p yields then 

PK(y) i-{i-pK{y)) J J 

that 



F(i- PK (y)) 

Pk(v) 



is bounded by (m-l)/3. Thus < 2(m-l)/3 

and we can appeal to bounded convergence as we take the limit in 
(4.17). With (4.18) and (4.19) we get 

2<m — 1 )B f K ° 

A(A) ~ c K 1 - )l \ sin 3 (7n//A )e^ dy, as K | A . 

Evaluating the integral gives 

12 (m-l)/3 7T 3 (e" A ° + l) 

Finally, A(A) ~ vr 2 (A - A )A7 3 as A | A which follows from the 
linearisation 

A(A) = (m-l)/3- /j, 2 /2 - n 2 /2K 2 

TV 2 7l 2 



2A 2 2A 2 



vr 2 (A-Ao + A ) 2 vr 2 A 2 



2A 2 A 2 2A 2 A 2 



7T 2 (A-A ) 7T 2 (A-A 



|2 



A A 2 2A 2 A 2 

and noting that the first term in the last line is the leading order term 
as A J, Kq. This completes the proof. □ 

5. Proof of Theorem 4. 

Proof of Theorem 4. Recall that (X, P D > K ) was defined as the 
process (X, P K ) conditioned on the event of survival and characterised 
via the change of measure in (3.7) and Theorem 5. 
Fix a A' > A and further denote by N t \^ K ^ the set of particles whose 
ancestors (including themselves) have not exited (0, A) up to time t. 
Then, for A < A', and for x G (0, A ) and A G 7 U 

lim P^? (A) = E« 1 A lim k °' K) -. 

KIKo x \ KIK Pk{x) 
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since Nt\, QK \ nas the same law under P K and P K . Suppose the par- 
ticles in ^Vtl(oftr) are ordered, for instance according to their spatial 
positions, and we write u\, ...,UN t \ . We can now expand the term 
within the expectation on the right-hand side as 

l-H ueNtkoK) (l-p K (x u (t))) 



Pk{x) 

<") = e ? ^#nd-f(^w)) 

i=l x j<i 

which is bounded from above by | N t \/ K \ \ {pk{x))~ x - Recall the asymp- 
totics for px in Theorem 3 and in particular Lemma 2, noting that these 
results hold uniformly in (0, K ). Since | N t \r QK \ | has finite expectation, 
we can apply the Dominated convergence theorem to the expression in 
(5.1), and we get 

PK{X Ut {t)) 



eKo ( 1a E sm{7TXu ^ t))/Ko)eflXuAt) 



sin(7rx/ ' Kq)c^ x 

Hence, for A G Jj, we arrive at 

hm e°> k {A) = e?° i A ^ ueNt . ; 0) 

KiK x y ' \ sm(7rx/K )e^ 

k ( Z Ko (t) 

where Z K ° is the martingale used in the change of measure in (2.2) 
in Section 2. The evolution under this change of measure is described 
in the paragraph following (2.2) and agrees with that of (X*, P*) as 
defined in Theorem 4. 

□ 



6. Super-Brownian motion in a strip. Recall from (1.1) that 

1j*1_ u A 

2dx 2 t^dxi 



the infinitesimal generator L is given as L = — /x-r-, x G (0, K), de 



fined for all functions u e C 2 (0, K) with u(0+) = u(K-) = 0. Chang- 
ing the domain to u E C 2 (0, K) with u"(0+) = u"(K-) = 0, then L 
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corresponds to Brownian motion with absorption (instead of killing) 
at and K. For technical reason, we will assume from now on that 
V K = {V^,t > 0} is the corresponding diffusion semi-group of Brow- 
nian motion with absorption and is therefore conservative. Note that 
all the results presented for branching Brownian motion with killing 
at and K also hold in the setting of absorption at and K when 
we restrict the process with absorption to particles within (0, K), in 
particular when defining N t as the number of particles who are alive 
and have not been absorbed at time t. 

Suppose Y = (Y t ,t > 0) is a Super-Brownian motion with associated 
semi-group V and branching mechanism ip of the form 

POO 

= —aX + /3X 2 + / (e~ Xy -l + Xy) Tl(dy), X > 0, 
Jo 

where a = —^'(0+) G (0, oo), (3 > and II is a measure concentrated 
on (0, oo) satisfying Jj Q ^(iAx 2 ) U(dx) < oo. For an initial configura- 
tion f] G Aif(0,K), the space of finite measures supported on (0,K), 
we denote the law of Y by Pff . The existence of this class of superpro- 
cesses follows from [Dyn91]. 

Since a = — ip'(0+) > 0, the function ip is the branching mecha- 
nism of a supercitcal continuous-state branching process (CSBP), say 
Z. We assume henceforth that ip satisfies the non-explosion condition 
Jo+ l^( s )l 1 ds = oo and further that ip(oo) = oo. The last condition, 
together with ?//(0+) < 0, ensures that ip has a unique positive root A*. 
The parameter A* is the survival rate of Z in the sense that the probabil- 
ity of the event {lim^oo Z t = 0} given Z = x is e~ x * x , which is strictly 
positive. We further assume from now on that f + °° (-^(s))" 1 ds < oo, 
which guarantees that the event {lim^oo Z t = 0} agrees with the event 
of extinction, that is {3t > : Z t = 0} a.s. This implies in turn that, for 
the Super-Brownian motion Y , the event of becoming extinguished and 
the event of extinction agree P -a.s. We denote the event of extinction 
of Y by 8 = {3t > : Y t (0, K) = 0}, where Y t (0, K) is the total mass 
within (0, K) at time t. We can characterise the P^ -superdiffusion via 
its Laplace functional (see e.g. Section 4.1.1 in [Dyn02]). 

Lemma 6. For all f G B+(0, K), 

E*(e~ if ' Yt) ) = e-^f^, 7] G M f (0,K), t > 0, 
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where uf(x,t) is the unique non-negative solution to the semi-group 
equation 

(6.1) Zf(x,t) = r*\f (•)](*)- fvl s [i>{uf{;s))]{x)ds. 

Jo 

We call the function u*f(x, t) the Laplace functional of (Y, P^). 

We have used the notation (/, 77) = f(x)r)(dx), for 77 G J\4f[0,K]. 
We define the survival rate wk of the P^-superdiffusion as the function 
satisfying 

Pf (£) = exp{-w K (x)}, for 77 6 M f [0,K), 

(6.2) 

and, taking / = 9 constant in Lemma 6 , we can deduce that 
- log P*{£) = lim hm(u e {-),ri) = (wk,v)- 

1 t-^-oo 8— too 

It can be derived, again by Lemma 6, that wk is a solution to 

(6.3) Lu - if>(u) = with u(0) = u(K) = 0. 

Analogous to Theorem 1, it is possible to give a necessary and sufficient 
condition for a positive survival rate which follows from a spine change 
of measure argument in the spirit of Section 2, now using the Pf- 
martingale 

(6.4) Z K (t) = / sin(7rx/K)e^-^ K)t Y t (dx), t > 0, 

Jo 

where here X(K) = —ip'(0+) — fi 2 /2 — tt 2 /2K 2 . Assuming henceforth in 
addition that x logxIl(dx) < 00, one can then show, in the fashion 
of Kyprianou et al. [KMSLR], that Z K is an L 1 (P^) martingale if 
and only if A (A') > 0. It can thus be concluded that wk is positive if 
X(K) > 0. 

Let us now establish the connection between the P^-superdiffusion and 
a P A -branching diffusion via the following relations. Set 

(6.5) F(s) = -^(A*(l-s)), s 6(0,1), 

(6.6) w K {x) = \*p K (x), xe(0,K), 
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where pk is the survival probability of the P -branching diffusion. 
Bertoin et al. [BFM08] show that (6.5) is the branching mechanism of 
a Galton- Watson process and they identify the Galton- Watson process 
with branching mechanism F of (6.5) as the backbone of the CSBP with 
branching mechanism ip. If we can show that wk in (6.6) is indeed the 
survival rate wk then the following Theorem is a direct consequence of 
Theorem 1 and Theorem 3. 

Theorem 7. (i) If fi < ^/-2-0'(O+) and K > K where K Q ■= 
7r(^/— 2ip'(0+))~ 1 , thenwK(x) > Oforallx G (0, K); otherwise wk(x) = 
for allx G [0,K\. 

(ii) Uniformly for x G (0, Kq), as K l Kq, 



Proof of Theorem 7. The relation in (6.5) gives (m - 1)(3 = 
F'(l—) = —if)'(0+) and hence the K in Theorem 7 is the same as the 
one in Theorem 1 and the X(K) defined earlier in this section agrees 
with X(K) as in Proposition 1. In particular, X(K) > if and only if 
fi < y/-2i/j'(0+) and K > K . 

Suppose [i < a/— 2ip'(0+) and K > K . By Remark 3, px is the unique 
non-trivial solution to Lu — F(l—u) = on (0, K) with u(0) = u(K) = 
0. Using (6.5) it follows then that wk given by (6.6) solves (6.3) and 
the uniqueness carries over. That is, for ji < \J — 2^'(0+) and K > K , 
wk is the unique non-trivial solution to (6.3). On the other hand, we 
know that wk solves (6.3) and, by the spine argument we mentioned 
after (6.4), we know that wk is positive within (0,K). By uniqueness, 
we have wk = wk- 

Suppose fi > \J— 2ip'{§+) or K < Kq. Then px is identically zero and 

(2.4) does not have a non-trivial solution. By the transformation in 

(6.5) , the same holds true for (6.3) and since wk is always a solution 
to (6.3) it must be equal to zero. Thus (6.6) holds true again. 

The result is now a direct consequence of Theorems 1 and 3. □ 

We will now outline the backbone decomposition for the P^-super 
-diffusion which consists of a copy of (Y, P^) conditioned on becom- 
ing extinct and further independent copies of (Y, P ) conditioned on 



w K (x) ~ X*{K-K ) 



(K V + 7r 2 )(K V + 97T 2 
12^'(0+)vrK3( e ^o + i) 



sm(iix/K )e' 



(6.7) 
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becoming extinct which immigrate along a P B '^-branching diffusion. 
Recall that the P B,it -branching diffusion is also the backbone of the 
P^-branching diffusion (Theorem 6). 

Let us begin by studying the process (Y, P K ) conditioned on becoming 
extinct. 

Proposition 8. Define for rj £ Mf[0, K) and t > 0, 



dP*> K 



dP* 



e -(w K ,Yt) 
e -(w K ,v) 

ft 



where (F t , t > 0) is the natural filtration generated by (Y,P^). Then 
(Y,P*' K ) is equal in law to (Y, P*(-\£ )). Further (Y,P^ K ) has spa- 
tially dependent branching mechanism 

4j r ' k {s, x) = 4j(s + w K (x)) - 4j(w k (x)), s > and x e [0, K] 

and diffusion semigroup V K . 

The proof of Proposition 8 is just a straightforward adaptation of the 
proof of Lemma 2 in [BKMS11] and thus omitted. We point out that 
the motion of the P^' - superdiffusion remains unchanged and it is 
therefore different from the motion of the P R) ^-branching diffusion in 
Definition 4. However, set uf(x,t) = A*(l —pK(x))(l — uf(x,t)), where 
is the Laplace functional of the P^'^-branching diffusion (Laplace 
functionals for branching diffusions are defined in a similar fashion 
to Lemma 6, see for instance Section 4.1.4 [Dyn02]). Then together 
with the relations (6.5) and (6.6) we can find that uf is the Laplace 
functional of the P R '^-superdiffusion. Thus the P R,i:f -superdiffusion 
can in this way be seen as the analogue of the P^'^-branching diffusion 
in the superdiffusion setting. 

We need to introduce some more notation before we can establish the 
backbone decomposition. Associated to the laws {Pf^ K ,x £ [0,-fT]} 
is the family of the so-called excursion measures {N R,K ,x £ [0,-fT]}, 
defined on the same measurable space, which satisfy 

(6.8) N£*(l - exp{-</, F t )}) = -logEf/(exp{-(/,y t )}), 

for any / £ B + [0, K] and t > 0. These measures are formally defined 
and studied in Dynkin and Kuznetov [DK04]. Intuitively speaking, the 
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branching property implies that P^ ,K is an infinitely divisible measure 
on the path space of Y and (6.8) is a 'Levy-Khinchine' formula in which 
N B,K plays the role of the Levy measure. In this sense, N B,K can be con- 
sidered as the rate at which P^'^-superdiffusions with infinitesimally 
small initial mass contribute to a unit mass at position x. Further we 
define, for n > 2, x £ (0, K), 



Pn{dy,x) 



(3w K (x) 2 5 (dy)l {n=2} + w K (x) n ^e w ^yU(dy) 
qn' K (x)w K (x)/3 B > K (x) 



which will turn out to be the distribution of the initial mass of the 
immigrating P^'^-superdiffusions at branch points of the backbone on 
the event of n offspring. 

Definition 1. LetK > K and v £ M a (0, K). LetX B = (X B ,t > 
0) be a P B ' K - branching diffusion with initial configuration v . Suppose 
I NR ' K = (lf' K ,t > 0), I pR ' K = {lf R ' K ,t > 0) and P = {I p t) t > 0) 
are three immigration processes (defined below) which are, condition- 
ally on X B , independent of each other. Then we define the process 
Y D = (Y t D ,t>0) by 

Y t D = lf' K + If*'* + II t > 
and denote its law by P^ ,K . 

The immigration processes are constructed as follows: 

(i) Continuous immigration: The process if ' = (if ' ,t > 0) is 
defined as 

r"=E E y ^ ,s \ t>o, 

ueT B r£ At<s<a B At 

where, given X B , independently for each u £ T B such that r B < t, the 
processes are countable in number and correspond to Poissonian 

immigration along the space-time trajectory {(x B (s),s) : s £ (T B ,a B ]} 
with rate 2Bds x dN R B , s . 

(ii) Discontinuous immigration: The process If ' = (7 t p ' ,t > 0) is 
defined as 

E Y ^ s \ t>o, 

u£T B t, b At<s<a B At 
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where, given X , independently for each u G T such that r R < t, the 
processes Y ( 2 >"' s ) are countable in number and correspond to Poissonian 
immigration along the space-time trajectory {(x B (s),s) : s G (T B ,a B ]} 
with rate ds x J °° yexp{-w K (x B {s))y} U(dy) x dPfy B . 

(Hi) Immigration at branch points: The process I p = > 0) is 

defined as 

where, given X B , independently for each u G T B stxc/i £/iat < i, 
Y (3.") ^ s an independent copy of (Y, Py ^ ) issued at space-time 

position (x B (a u ),a u ). At a branch point of u with n > 2 offspring the 
initial mass Y u is distributed according to p n (dy,x B (o~ u )). 

Theorem 8 (Backbone decomposition). For K > K and 77 G 
M f [0,K], let Y R = (Y t R ,t > 0) be an independent copy of (Y,P R < K ). 
Suppose that v is a Poisson random measure on (0, K) with intensity 
WK{x)rj(dx). Let (Y D ,P® ,K ) be the process constructed in Definition 
1. Define the process Y = (Y t ,t > 0) by 

(6.9) Y t = Y* + Y t D , t>0, 

and denote its law by P^ . Then the process (Y, P^) is Markovian and 
equal in law to (Y,P^). 

The proof of Theorem 8 is a simple adaptation of the proofs of Theo- 
rem 1 and 2 in [BKMS11] and therefore omitted. In principle it should 
be possible to give a proof analogous to the one presented in Section 
3, using that (Y, Pff) conditioned on non-extinction arises from the 
martingale change of measure 

dP u K 

and showing that (Y, Q^) agrees in law with the process (Y D , P®' K ) 
of Definition 1. 

The analogy between the P^-branching diffusion and the P^-super- 
diffusion indicates that there is a quasi- stationary limit result equiva- 
lent to Theorem 4. Let us begin with constructing the analogue of the 



1 _ e ~{wK,Yt) 
1 _ p -{w K ,u) 



t > 0, 
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process (X*, P*) in Theorem 4 for the superdiffusion setting. We intro- 
duce the family of N- measures now associated with the laws (P^, l \ x G 
[0, K Q ]). Consider the family {Nf°,x G [0, K ]} satisfying 

N?>(1 - exp (/, Y t )) — — log £f(e-<^>), 

for / G B+(0,K), t > 0. 

Let r) G .M f (0, K). Suppose £* = (£*, t > 0) is a Brownian motion con- 
ditioned to stay in (0, K ) with initial position x distributed according 
to 

(e.io) si . n ^.y; , e( o,if„). 

Let /f Ko = > 0) and i^ Ko = (lf K \t > 0) be two immigration 

processes (defined below) which, conditionally on £*, are independent 
of each other. Then we define the process Y s = (Y t s ,t > 0) by 

(6.11) Y t s = lf + lf K \ t>0. 



The immigration processes I nK ° and I pK ° are defined pathwise as fol- 
lows. 

(i) Continuous immigration: The process If = (If , t > 0) is defined 

as 



where, given £*, the processes F (N ' s) are countable in number and 
correspond to Poissonian immigration along the space-time trajectory 
{(€> s) : s > 0} with rate 2/3c/s x 

(ii) Discontinuous immigration: The process J pA ° = (lf K ° , £ > 0) is 
defined as 



s<t 



where, given £*, the processes y( pK °' s ) are countable in number and 
correspond to Poissonian immigration along the space-time trajectory 
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{(£, s) : a > 0} with ds x / °° y n(dy) x . 
Then define the process Y* = (Y t *,t > 0) by setting 

(6.12) Y*=Yl + Yf, t>0, 

where Y' is an independent copy of (Y, P^ )- We denote the law of Y* 
by P*. The evolution of Y* under P* can thus be seen as a path- wise 
description of Evans' immortal particle picture in [Eva93] for the crit- 
ical width Kq] for a similar construction of Evans' immortal particle 
picture see Kyprianou et al. [KMSLR]. 

Further, we note that (Y*, P^°) has the same law us Y under the mea- 
sure which has martingale density Z K °(t) of (6.4) with respect to P^°; 
for similar results see for instance Englander and Kyprianou [EK04], 
Kyprianou et al. [KMSLR] and Liu et al. [LRS09]. 

Theorem 9. Let K > K and 77 e Mf[0,K ]. For a fixed time 
t > 0, the law ofY t under the measure limiqjc P„{- \ lim^oo \ \Y t \ \ > 0) 
is equal to Y t * under P* . 

Sketch of the proof of Theorem 9. By Theorem 8, (Y,P^) 
is equal in law to (Y, P5). The latter is equal in law to (Y, Q^) where 



1 _ e -(w K ,r)) 
Ft 



(w K ,Yt) 

t > 0. 



The uniform asymptotics for wk in Theorem 7 let us conclude that 

l_ e -(^> (w K ,Y t ) 
urn = Iim 

KIK \- e -{w K ,ri) KiKo (w K ,v} 

f*° smi-Kx/Koje^ Y t {dx) _ Z Ka {t) 
J Ko sm(7rx/K )e^ fi(dx)) ~ P^O)' 

where Z Ko is the martingale in (6.4). As mentioned before, the law of 
Y under a change of measure with Z Ko is equal to (Y*, P^°). □ 
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